Robustness of cancer microbiome signals over a broad range of methodological variation.

In 2020, we identified cancer-specific microbial signals in The Cancer Genome Atlas (TCGA) [1]. Multiple peer-reviewed papers independently verified or extended our findings [2-12]. Given this impact, we carefully considered concerns by Gihawi et al. [13] that batch correction and database contamination with host sequences artificially created the appearance of cancer type-specific microbiomes. (1) We tested batch correction by comparing raw and Voom-SNM-corrected data per-batch, finding predictive equivalence and significantly similar features. We found consistent results with a modern microbiome-specific method (ConQuR [14]), and when restricting to taxa found in an independent, highly-decontaminated cohort. (2) Using Conterminator [15], we found low levels of human contamination in our original databases (~1% of genomes). We demonstrated that the increased detection of human reads in Gihawi et al. [13] was due to using a newer human genome reference. (3) We developed Exhaustive, a method twice as sensitive as Conterminator, to clean RefSeq. We comprehensively host-deplete TCGA with many human (pan)genome references. We repeated all analyses with this and the Gihawi et al. [13] pipeline, and found cancer type-specific microbiomes. These extensive re-analyses and updated methods validate our original conclusion that cancer type-specific microbial signatures exist in TCGA, and show they are robust to methodology.

Signatures of HIV and Major Depressive Disorder in the Plasma Microbiome.

Inter-individual differences in the gut microbiome are linked to alterations in inflammation and blood-brain barrier permeability, which may increase the risk of depression in people with HIV (PWH). The microbiome profile of blood, which is considered by many to be typically sterile, remains largely unexplored. We aimed to characterize the blood plasma microbiome composition and assess its association with major depressive disorder (MDD) in PWH and people without HIV (PWoH). In this cross-sectional, observational cohort, we used shallow-shotgun metagenomic sequencing to characterize the plasma microbiome of 151 participants (84 PWH and 67 PWoH), all of whom underwent a comprehensive neuropsychiatric assessment. The microbial composition did not differ between PWH and PWoH or between participants with MDD and those without it. Using the songbird model, we computed the log ratio of the highest and lowest 30% of the ranked classes associated with HIV and MDD. We found that HIV infection and lifetime MDD were enriched in a set of differentially abundant inflammatory classes, such as Flavobacteria and Nitrospira. Our results suggest that the circulating plasma microbiome may increase the risk of MDD related to dysbiosis-induced inflammation in PWH. If confirmed, these findings may indicate new biological mechanisms that could be targeted to improve treatment of MDD in PWH.

Oral mucosal breaks trigger anti-citrullinated bacterial and human protein antibody responses in rheumatoid arthritis.

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

Standardized multi-omics of Earth's microbiomes reveals microbial and metabolite diversity.

Despite advances in sequencing, lack of standardization makes comparisons across studies challenging and hampers insights into the structure and function of microbial communities across multiple habitats on a planetary scale. Here we present a multi-omics analysis of a diverse set of 880 microbial community samples collected for the Earth Microbiome Project. We include amplicon (16S, 18S, ITS) and shotgun metagenomic sequence data, and untargeted metabolomics data (liquid chromatography-tandem mass spectrometry and gas chromatography mass spectrometry). We used standardized protocols and analytical methods to characterize microbial communities, focusing on relationships and co-occurrences of microbially related metabolites and microbial taxa across environments, thus allowing us to explore diversity at extraordinary scale. In addition to a reference database for metagenomic and metabolomic data, we provide a framework for incorporating additional studies, enabling the expansion of existing knowledge in the form of an evolving community resource. We demonstrate the utility of this database by testing the hypothesis that every microbe and metabolite is everywhere but the environment selects. Our results show that metabolite diversity exhibits turnover and nestedness related to both microbial communities and the environment, whereas the relative abundances of microbially related metabolites vary and co-occur with specific microbial consortia in a habitat-specific manner. We additionally show the power of certain chemistry, in particular terpenoids, in distinguishing Earth's environments (for example, terrestrial plant surfaces and soils, freshwater and marine animal stool), as well as that of certain microbes including Conexibacter woesei (terrestrial soils), Haloquadratum walsbyi (marine deposits) and Pantoea dispersa (terrestrial plant detritus). This Resource provides insight into the taxa and metabolites within microbial communities from diverse habitats across Earth, informing both microbial and chemical ecology, and provides a foundation and methods for multi-omics microbiome studies of hosts and the environment.

Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Multiomic Analyses of Nascent Preterm Infant Microbiomes Differentiation Suggest Opportunities for Targeted Intervention.

The first week after birth is a critical time for the establishment of microbial communities for infants. Preterm infants face unique environmental impacts on their newly acquired microbiomes, including increased incidence of cesarean section delivery and exposure to antibiotics as well as delayed enteral feeding and reduced human interaction during their intensive care unit stay. Using contextualized paired metabolomics and 16S sequencing data, the development of the gut, skin, and oral microbiomes of infants is profiled daily for the first week after birth, and it is found that the skin microbiome appears robust to early life perturbation, while direct exposure of infants to antibiotics, rather than presumed maternal transmission, delays microbiome development and prevents the early differentiation based on body site regardless of delivery mode. Metabolomic analyses identify the development of all gut metabolomes of preterm infants toward full-term infant profiles, but a significant increase of primary bile acid metabolism only in the non-antibiotic treated vaginally birthed late preterm infants. This study provides a framework for future multi-omic, multibody site analyses on these high-risk preterm infant populations and suggests opportunities for monitoring and intervention, with infant antibiotic exposure as the primary driver of delays in microbiome development.

A comparison of six DNA extraction protocols for 16S, ITS and shotgun metagenomic sequencing of microbial communities.

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

The Host-Microbiome Response to Hyperbaric Oxygen Therapy in Ulcerative Colitis Patients.

BACKGROUND & AIMS: Hyperbaric oxygen therapy (HBOT) is a promising treatment for moderate-to-severe ulcerative colitis. However, our current understanding of the host and microbial response to HBOT remains unclear. This study examined the molecular mechanisms underpinning HBOT using a multi-omic strategy. METHODS: Pre- and post-intervention mucosal biopsies, tissue, and fecal samples were collected from HBOT phase 2 clinical trials. Biopsies and fecal samples were subjected to shotgun metaproteomics, metabolomics, 16s rRNA sequencing, and metagenomics. Tissue was subjected to bulk RNA sequencing and digital spatial profiling (DSP) for single-cell RNA and protein analysis, and immunohistochemistry was performed. Fecal samples were also used for colonization experiments in IL10-/- germ-free UC mouse models. RESULTS: Proteomics identified negative associations between HBOT response and neutrophil azurophilic granule abundance. DSP identified an HBOT-specific reduction of neutrophil STAT3, which was confirmed by immunohistochemistry. HBOT decreased microbial diversity with a proportional increase in Firmicutes and a secondary bile acid lithocholic acid. A major source of the reduction in diversity was the loss of mucus-adherent taxa, resulting in increased MUC2 levels post-HBOT. Targeted database searching revealed strain-level associations between Akkermansia muciniphila and HBOT response status. Colonization of IL10-/- with stool obtained from HBOT responders resulted in lower colitis activity compared with non-responders, with no differences in STAT3 expression, suggesting complementary but independent host and microbial responses. CONCLUSIONS: HBOT reduces host neutrophil STAT3 and azurophilic granule activity in UC patients and changes in microbial composition and metabolism in ways that improve colitis activity. Intestinal microbiota, especially strain level variations in A muciniphila, may contribute to HBOT non-response.

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy.

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

Salivary bacterial signatures in depression-obesity comorbidity are associated with neurotransmitters and neuroactive dipeptides.

BACKGROUND: Depression and obesity are highly prevalent, often co-occurring conditions marked by inflammation. Microbiome perturbations are implicated in obesity-inflammation-depression interrelationships, but how the microbiome mechanistically contributes to pathology remains unclear. Metabolomic investigations into microbial neuroactive metabolites may offer mechanistic insights into host-microbe interactions. Using 16S sequencing and untargeted mass spectrometry of saliva, and blood monocyte inflammation regulation assays, we identified key microbes, metabolites and host inflammation in association with depressive symptomatology, obesity, and depressive symptomatology-obesity comorbidity. RESULTS: Gram-negative bacteria with inflammation potential were enriched relative to Gram-positive bacteria in comorbid obesity-depression, supporting the inflammation-oral microbiome link in obesity-depression interrelationships. Oral microbiome was more highly predictive of depressive symptomatology-obesity co-occurrences than of obesity or depressive symptomatology independently, suggesting specific microbial signatures associated with obesity-depression co-occurrences. Mass spectrometry analysis revealed significant changes in levels of signaling molecules of microbiota, microbial or dietary derived signaling peptides and aromatic amino acids among depressive symptomatology, obesity and comorbid obesity-depression. Furthermore, integration of the microbiome and metabolomics data revealed that key oral microbes, many previously shown to have neuroactive potential, co-occurred with potential neuropeptides and biosynthetic precursors of the neurotransmitters dopamine, epinephrine and serotonin. CONCLUSIONS: Together, our findings offer novel insights into oral microbial-brain connection and potential neuroactive metabolites involved.

Swapping Metagenomics Preprocessing Pipeline Components Offers Speed and Sensitivity Increases.

Increasing data volumes on high-throughput sequencing instruments such as the NovaSeq 6000 leads to long computational bottlenecks for common metagenomics data preprocessing tasks such as adaptor and primer trimming and host removal. Here, we test whether faster recently developed computational tools (Fastp and Minimap2) can replace widely used choices (Atropos and Bowtie2), obtaining dramatic accelerations with additional sensitivity and minimal loss of specificity for these tasks. Furthermore, the taxonomic tables resulting from downstream processing provide biologically comparable results. However, we demonstrate that for taxonomic assignment, Bowtie2's specificity is still required. We suggest that periodic reevaluation of pipeline components, together with improvements to standardized APIs to chain them together, will greatly enhance the efficiency of common bioinformatics tasks while also facilitating incorporation of further optimized steps running on GPUs, FPGAs, or other architectures. We also note that a detailed exploration of available algorithms and pipeline components is an important step that should be taken before optimization of less efficient algorithms on advanced or nonstandard hardware. IMPORTANCE In shotgun metagenomics studies that seek to relate changes in microbial DNA across samples, processing the data on a computer often takes longer than obtaining the data from the sequencing instrument. Recently developed software packages that perform individual steps in the pipeline of data processing in principle offer speed advantages, but in practice they may contain pitfalls that prevent their use, for example, they may make approximations that introduce unacceptable errors in the data. Here, we show that differences in choices of these components can speed up overall data processing by 5-fold or more on the same hardware while maintaining a high degree of correctness, greatly reducing the time taken to interpret results. This is an important step for using the data in clinical settings, where the time taken to obtain the results may be critical for guiding treatment.

The microbiome and prostate cancer.

There is growing evidence that the microbiome is involved in development and treatment of many human diseases, including prostate cancer. There are several potential pathways for microbiome-based mechanisms for the development of prostate cancer: direct impacts of microbes or microbial products in the prostate or the urine, and indirect impacts from microbes or microbial products in the gastrointestinal tract. Unique microbial signatures have been identified within the stool, oral cavity, tissue, urine, and blood of prostate cancer patients, but studies vary in their findings. Recent studies describe potential diagnostic and therapeutic applications of the microbiome, but further clinical investigation is needed. In this review, we explore the existing literature on the discovery of the human microbiome and its relationship to prostate cancer.

Utilizing stability criteria in choosing feature selection methods yields reproducible results in microbiome data.

Feature selection is indispensable in microbiome data analysis, but it can be particularly challenging as microbiome data sets are high dimensional, underdetermined, sparse and compositional. Great efforts have recently been made on developing new methods for feature selection that handle the above data characteristics, but almost all methods were evaluated based on performance of model predictions. However, little attention has been paid to address a fundamental question: how appropriate are those evaluation criteria? Most feature selection methods often control the model fit, but the ability to identify meaningful subsets of features cannot be evaluated simply based on the prediction accuracy. If tiny changes to the data would lead to large changes in the chosen feature subset, then many selected features are likely to be a data artifact rather than real biological signal. This crucial need of identifying relevant and reproducible features motivated the reproducibility evaluation criterion such as Stability, which quantifies how robust a method is to perturbations in the data. In our paper, we compare the performance of popular model prediction metrics (MSE or AUC) with proposed reproducibility criterion Stability in evaluating four widely used feature selection methods in both simulations and experimental microbiome applications with continuous or binary outcomes. We conclude that Stability is a preferred feature selection criterion over model prediction metrics because it better quantifies the reproducibility of the feature selection method.

Class Ib MHC-Mediated Immune Interactions Play a Critical Role in Maintaining Mucosal Homeostasis in the Mammalian Large Intestine.

Lymphocytes within the intestinal epithelial layer (IEL) in mammals have unique composition compared with their counterparts in the lamina propria. Little is known about the role of some of the key colonic IEL subsets, such as TCRαß+CD8+ T cells, in inflammation. We have recently described liver-enriched innate-like TCRαß+CD8αα regulatory T cells, partly controlled by the non-classical MHC molecule, Qa-1b, that upon adoptive transfer protect from T cell-induced colitis. In this study, we found that TCRαß+CD8αα T cells are reduced among the colonic IEL during inflammation, and that their activation with an agonistic peptide leads to significant Qa-1b-dependent protection in an acute model of colitis. Cellular expression of Qa-1b during inflammation and corresponding dependency in peptide-mediated protection suggest that Batf3-dependent CD103+CD11b- type 1 conventional dendritic cells control the protective function of TCRαß+CD8αα T cells in the colonic epithelium. In the colitis model, expression of the potential barrier-protective gene, Muc2, is enhanced upon administration of a Qa-1b agonistic peptide. Notably, in steady state, the mucin metabolizing Akkermansia muciniphila was found in significantly lower abundance amid a dramatic change in overall microbiome and metabolome, increased IL-6 in explant culture, and enhanced sensitivity to dextran sulfate sodium in Qa-1b deficiency. Finally, in patients with inflammatory bowel disease, we found upregulation of HLA-E, a Qa-1b analog with inflammation and biologic non-response, in silico, suggesting the importance of this regulatory mechanism across species.

Efficient computation of Faith's phylogenetic diversity with applications in characterizing microbiomes.

The number of publicly available microbiome samples is continually growing. As data set size increases, bottlenecks arise in standard analytical pipelines. Faith's phylogenetic diversity (Faith's PD) is a highly utilized phylogenetic alpha diversity metric that has thus far failed to effectively scale to trees with millions of vertices. Stacked Faith's phylogenetic diversity (SFPhD) enables calculation of this widely adopted diversity metric at a much larger scale by implementing a computationally efficient algorithm. The algorithm reduces the amount of computational resources required, resulting in more accessible software with a reduced carbon footprint, as compared to previous approaches. The new algorithm produces identical results to the previous method. We further demonstrate that the phylogenetic aspect of Faith's PD provides increased power in detecting diversity differences between younger and older populations in the FINRISK study's metagenomic data.

EMPress Enables Tree-Guided, Interactive, and Exploratory Analyses of Multi-omic Data Sets.

Standard workflows for analyzing microbiomes often include the creation and curation of phylogenetic trees. Here we present EMPress, an interactive web tool for visualizing trees in the context of microbiome, metabolome, and other community data scalable to trees with well over 500,000 nodes. EMPress provides novel functionality-including ordination integration and animations-alongside many standard tree visualization features and thus simplifies exploratory analyses of many forms of 'omic data.IMPORTANCE Phylogenetic trees are integral data structures for the analysis of microbial communities. Recent work has also shown the utility of trees constructed from certain metabolomic data sets, further highlighting their importance in microbiome research. The ever-growing scale of modern microbiome surveys has led to numerous challenges in visualizing these data. In this paper we used five diverse data sets to showcase the versatility and scalability of EMPress, an interactive web visualization tool. EMPress addresses the growing need for exploratory analysis tools that can accommodate large, complex multi-omic data sets.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment.

BACKGROUND: Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical-grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients including COVID-19+ patients. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent. RESULTS: Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4× higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic (SYN) swab. The limit of detection (LoD) of SARS-CoV-2 from floor samples collected using the consumer-grade plastic (CGp) or research-grade plastic The Microsetta Initiative (TMI) swabs was similar or better than the SYN swab, further suggesting that swab type does not impact RNA recovery as measured by the abundance of SARS-CoV-2. The LoD for TMI was between 0 and 362.5 viral particles, while SYN and CGp were both between 725 and 1450 particles. Lastly microbiome analyses (16S rRNA gene sequencing) of paired samples (nasal and floor from same patient room) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type, but instead driven by the patient and sample type. CONCLUSIONS: Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer-grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity and lack of antibiotics in these samples makes it possible to perform concomitant microbiome analyses. Video abstract.

Potential Immune-Microbiome Interactions in Breast Cancer May Advance Treatment: What's Holding Us Back?

The impact of the human microbiome, the diverse collection of microorganisms that inhabit nearly every niche in the human body, in shaping the immune response to dysbiotic events is apparent if poorly understood, particularly in complex, evolving disease states such as breast cancer. The impacts can be both indirect via metabolites and immune-interactions along the skin, gut, and oral cavities where the microbial communities are most abundant, or direct in the tumor microenvironment where microbial activities can promote growth or clearance of cancerous cells. Based on reports of using gut microbial signatures to predict therapeutic efficacy, the role that gut microbes and their metabolites may play in shaping the success or failure of immunotherapy has been extensively reviewed. In this review, we dissect the evidence for the direct and distal impact of microbes on oncogenesis, tumor growth and the immune responses to combat or promote tolerance of breast cancer tumors. Implementation of robust, valid analyses and methods are lacking in the field, and we provide recommendations for researchers and clinicians to work together to characterize the micro-biome-immune-breast cancer interactions that will hopefully enable the next generation of biomarkers and targets for improving disease outcomes.